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Abstract. We show non-integrability of the nonhnear lattice of Fermi-Pasta-Ulam type 
via singularity analysis of normal variational equations of Lame type. 



1. From a Nonlinear Lattice to Lame Equations 

We consider the following one-dimensional lattice: 

n -j^ n+1 

^ i=i ^ i=i 

where 

^(X) = ^X2 + ^X4 + ... + ^X2-. (2) 
2 4 2m 

Fermi-Pasta-Ulam(FPU) lattice [2] is a special type of the systems with the potential 
function (2) as follows: 

n n+1 n+1 

i=l i=l i=l 

If we impose the fixed boundary condition as 

qo = qn+i =0, n = odd, (4) 

it is easy to check that 

n — 1 

T: qi = C^{t),q2=0,q3 = -CcP{t),---,qn-i=0,qn = {-l) — CcP{t) (5) 

is a special solution. Thus, the equation of (p{t) is equivalent to the following Hamiltonian 
system with one degree of freedom: 
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where Hamiltonian is 

H{cf>, 4>) = \{4>? + ^^24>'' + ^<t>' + ■■■ + = Const. (7) 

Then the total energy e is given by 

6 = = <^)^C2 = ^C2(/X2 + ^M4C2 + . . . (8) 
for the initial condition (5). In the case of the FPU lattice, we can determine C as follows: 



C = 



By combining (7) with (8), the underlying equation of (f){t) can be rewritten by the differ- 
ential equation of <p{t) as 

= 72(1 - + |(i _ ^4) + . . . + _ ^2n.^^ (10) 

where 

72m(e, {/^2, |j = 1, • • • , m}) = /X2mC2"^"'- (H) 

In the case of the FPU lattices (3), the solution of this differential equation (10) with 
the condition 

72m=4 / (12) 

is given explicitly by the elliptic function 

^{t) =cn{k;at), (13) 

where 

a = V272 + 274, k = J , (14) 

V 272 + 274 

cn{k; at) is the Jacobi cn elliptic function, and k is the modulus of the elliptic integral. 
We remark that because 



72 + 74 = /^2 + C'^A«4 = -iM + — -tM4 > 0, (15) 

y n + 1 

holds for //4 > 0, jLt2 > 0, the modulus of the elliptic function k satisfies the following 
relation: 

< < (16) 

Thus, the special solutions of the FPU lattices for )U4 > 0, /X2 > have the two fundamental 
periods in the complex time plane as follows: 

= r.(,.) = ?m±^, (17) 

a a 
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where K(k) and K'{k) are the complete elhptic integrals of the first kind: 



K{k) 



dv 



^/{l-v''){l-k'^v'^ 



= , K'ik) 



dv 



^{i-v'){i-{i-k'')v''y 



(18) 



Poles are located at t = r, where r = ^^^^ + i^-^ (mod Ti,T2) in the parallelogram 
of each period cell. Let us consider the variational equations along these special solutions. 
The variational equations are obtained by 



^3 =tj = -J2 



k=l dqkdqj 



-(72 + + 5^6^ + • • • + (2m - l)72„<^2m-2)(2^^. - - ^.^^) for 1 < j < n, 



where = Cn+i = Vo = Vn+i = and = Sqj,r]j = dpj (1 < j < n). 
Moreover, these linear variational equations in the form of the vector 



d^ 
dt^^ 



(19) 
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(20) 



can be decoupled as follows. After we note that the eigenvalues of the n x n symmetric 
matrix 
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(21) 



are obtained as {4sin^(2^^^jy)|l < J < n.} by a normal orthogonal transformation G 
OGO~^, the variational equations (19) are rewritten in the decoupled form: 



-)(72 + 374<^' + • • • + (2m - l)72m</'""-')e;-(i) (1 < i < n), (22) 



^2(n + l)' 

where $f = O^. Clearly, these equations are written in the form of vector Hill's equation[3\ 



df^ 



+ A{t)i' = Q, A{t + T) = A{t), 



(23) 



where T = Ti,T2 in the case of m = 2. For j = we have the relation 



^ (ei-6 + ^5 + --- + (-l)'^en)- 



/ n+ 1 

Thus, the corresponding variational equation 

i'n+i = -2(72 + 37402 + • • • + (2m - l)72m0"""')CUi («) 



(24) 



(25) 
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has a time-dependent integral I{^,^;t) = I{^,r];t) because 
I{^,rj;t) = DH={iT^+^-^)H = ri-p + ^-Vg 

= C^{rii - 773 + % + ••• + i-l)'^Vn) (26) 

+2(C72</' + C74^' + • • • + C72m0'"^-')(a - 6 + ^5 + • • • + (-l)'^en), 



where 



= <^(^i-^3 + --- + (-i)^e;) 

+2</.(72 + 374</>2 + • • • + (2m - l)72m<^'"-')(ei - 6 + • • • + (-l)'^en) = 0. 



(27) 



We call Eq. (25) the tangential variational equation. On the other hands, a (2n — 2)- 
dimensional normal variational equation{NYE) is given by the equation of (22) with the 
tangential variational equation (25) removed as follows: 

V'j = -4sm2(2(i^)(72 + 374</'' + • • • + (2m - l)72m</'''"-')ei, .28) 
i'j = r]'^ for 1 < j(7^ ^) < n. ^ ' 

In case of the FPU lattice, the normal variational equation (28) becomes the Lame equation 
[7] 

^ - {Eisn''{k; at) + E2)y = 0, (29) 
where Ei = 12^^sin^(2^^^jy) and E2 are constants. 

2. Non-integrability Theorem 

Morales and Simo obtained the following theorem on the non-integrability based on the 
application of Picard-Vessiot theory to Ziglin's analysis[9, 10] for Hamiltonian systems 
with two degrees of freedom. 

Theorem 1 (Morales and Simo [4], 1994) When the normal reduced variational equa- 
tion is of Lame type, if A = E^a'^k'^ ^ m{m-\-l),m G N and the Lame equation satisfying 
this condition on A is not algebraically solvable(Brioschi-Halphen-Crawford and Baldas- 
sarri solutions), then the the initial Hamiltonian system does not have a first integral, 
meromorphic in a connected neighborhood of the integral curve T, which is functionally 
independent together with H. 

In case of the present analysis, A is given by the following formula: 

A = E,a'k' = 12sin2(^^^) = 6(1 - cos(^)). (30) 

We can easily check that cos(^) ^ Q if and only if j ^ ^^^}. When A(^Q, 

the above condition on the algebraic solvability of the Lame equation is not satisfied. Thus, 
to check the algebraic solvability of the Lame equations 
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it is sufficient to examine the following two cases: 

A = 6(l-cos(^)) = 3, A = 6(l-cos(— )) =9. (32) 

It is known [1] that the condition on A for the Brioschi-Halphen-Crawford solutions is 
given by 

A = m{m + 1), m + ^GAT, (33) 
and that the condition on A for the Baldassarri solutions is given by 

A = m(m + 1), m+-£-Z(J^ZU-Z\Z. (34) 
^ " 2 3 4 5 

However, the following relations 

m{m + 1) = 3 ^ m = ^ Q, m{m + 1) = 9 ^ m = ^ Q (35) 

hold, which guarantee that all n — 1 Lame equations (31) do not belong to the solvable 
case. In case of the systems with n degrees of freedom, we have ra — 1 Lame equations 
which corresponds to n — 1 normal variational equations. 

Thus, according to the steps in Ref. [4] we obtain the following theorem: 

Theorem 2 The FPU lattice for ^4 > 0, //2 > does not have n — 1 first integrals, 
meromorphic in a connected neighbourhood of the integral curve T, which are functionally 
independent together with H. 

We remark here that this theorem on the non-integrability does not depend on the total 
energy in contrast with the result about the non-integrability proof of the FPU lattice in 
the low energy limit [6] based on non-resonance checking [5] and the result about the non- 
integrability of the FPU lattice in the high energy limit based on the Kowalevski exponents 
of the homogeneous systems [8]. Here, it is conjectured that more general nonhomogeneous 
nonlinear lattice (1) would be also non-integrable in the sense of the present analysis. 
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